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ABSTRACT 

The  analytical  solution  of  the  heat  conduction  equation  for  a 
gun  tube  (circular,  thick-walled  cylinder)  assumed,  for  simplicity, 
to  be  of  infinite  length,  is  obtained.  The  temperature  function  is 
assumed  to  have  axial  and  angular  symmetry  and,  therefore,  is  depen¬ 
dent  only  upon  time  and  radial  distance.  The  axial  flow  has  been 
neglected  because  the  projectile  moves  through  the  tube  at  a  much 
faster  rate  than  heat  can  penetrate  the  tube.  The  initial  and 
boundary  conditions  have  been  considered  a?  close  as  possible  to 
those  encountered  during  the  repeated  firing  of  a  gun.  The  tempera¬ 
ture  distribution  in  a  barrel  can  be  determined  provided  that  the 
propellant  gas  temperature  during  the  firing  of  each  round  is  known 
and  that  the  nature  of  the  heat  transfer  from  the  gas  to  the  barrel 
wall  is  known. 
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NOMENCLATURE 


a  «  inside  radius  of  gun  tube 
b  =  outside  radius  of  gun  tube 
t  *  time 

hj,h2  =  gas  and  natural  convection  coefficients  of  gun  tube, 
respectively 

K  =  thermal  conductivity  of  gun  tube 

r  =  radial  distance  of  a  point  from  the  axis  of  the  cylinder 

T(r,t)  -  temperature  of  point  r  at  time  t 

Tg(t)  =  propellant  gas  temperature  function 

Ta(r,t)  =  a  transient  part  of  temperature  of  the  total  temperature 

Tb(r,t)  =  n  transient  part  of  temperature  of  the  total  temperature 

Ts(r)  =  a  steady-state  part  of  temperature  of  the  total  temperature 

Ty\  =  ambient  temperature 

a  =  thermal  diffusivity  of  barrel  material 

Bn  '  eigenvalues  of  equation  (30) 

J0,Y0  =  Bessel  functions  of  first  and  second  kinds  of  order  zero, 
respectively 

J,  Yx  =  Bessel  functions  of  first  and  second  kinds  of  order  one, 
respectively 

!0,K0  =  modified  Bessel  functions  of  first  and  second  kinds  of 
order  zero,  respectively 

I.jKj  =  modified  Bessel  functions  of  first  and  second  kinds  of 
order  one,  respectively 

f(r)  =  initial  temperature  distribution 


A 


INTRODUCTION 


The  solution  to  the  initial -boundary  value  problem  for  the 
conduction  of  heat  in  a  thick-walled  cylinder  is  of  considerable 
technological  importance,  particularly,  in  the  gun  tube  heat 
transfer  analysis.  These  problems  have  consequently  attracted 
considerable  attention,  and  a  number  of  special  solutions1-6 
have  been  developed.  M,ow,  Pascual,  and  Pascale1  dealt  with  the 
problem  of  radial  flow  of  heat  in  a  gun  tube  assumed  to  be 
infinitely  long.  In  their  investigations,  the  initial  and  boundary 
conditions  were  assumed  to  be  T(r,t)  =  0  for  t  -  0,  T(a,t)  =  <j>(t) 
for  t  *  0,  and  T(b,t)  =  0  for  all  t.  Apparently,  the  surface 
conductance  on  both  inner  and  outer  surfaces  of  a  gun  tube  has 
been  disregarded  in  their  analysis.  A  similar  problem,  with  a 
convective  boundary  condition  on  the  outer  surface  of  a  gun  tube, 
was  solved  by  Pascual,  Zv/eig,  and  Sutherland.2  By  virtue  of  the 
fact  that  the  initial  condition  was  assumed  to  be  zero  in  the  above 
two  investigations,  the  problem  solved  by  these  investigators  is 
essentially  the  first  round  problem.  Comenetz3  gave  a  solution  for 
the  temperature  distribution  in  a  hollow  cylinder  of  infinite  length 
supplied  with  heat  through  the  inner  surface.  The  rate  of  heat 
input  was  permitted  to  vary  linearly  (or  at  most  quadratically)  with 
time  and  the  initial  temperature  of  the  hollow  cylinder  was  assumed 
to  be  zero.  A  more  general  solution  of  heat  flowjn  a  finite,  hollow, 
circular  cylinder  was  given  by  Olcer.4  However,  Olcer's  method  is 
very  time-consuming  and  difficult  to  obtain  desirable  numerical  results, 
fu  solve  the  gun  barrel  heat-transfer  problem  practically,  a  relatively 
simple  solution  that  closely  approximates  the  actual  firing  conditions 
of  the  gun  is  necessary.  This  report  is  intended  to  serve  this 
purpose. 

OBJECTIVE 


One  of  the  many  causes  of  gun  barrel  failure  occurs  when  the 
effective  thermal  stress  level  in  the  interior  of  the  small  caliber 
gun  tube  exceeds  the  yield  strength  of  the  barrel  material. 

To  determine  the  thermal  stresses  in  a  barrel,  the  temperature 
distribution  within  a  barrel  must  first  be  determined.  Therefore, 
the  purpose  of  this  investigation  is  to  obtain  the  analytical  solution 
of  the  heat  conduction  equation  for  a  gun  tube  subjected  to  an  arbi¬ 
trary  initial  condition  and  to  boundary  conditions  that  are  as  close 
as  possible  to  those  encountered  during  the  actual  single  and  repeated 
fi ring. 


FORMULATION  OF  THF  PROBLEM 


The  unsteady  temperature  field  in  a  homogeneous,  isotropic  gun 
tube  is  assumed  to  have  axial  and  angular  symmetry  and,  therefore, 
is  dependent  only  upon  time  and  radial  distance.  For  simplicity,  the 
gun  tube  is  assumed  to  be  of  infinite  length,  and  the  thermal  prop¬ 
erties  assumed  to  be  independent  of  the  temperature  for  each  small 
time-interval  during  each  firing  round.  Then  the  heat  flow  in  a  gun 
tube  is  governed  by  the  well-known  Fourier  heat  conduction  equation.5 


a?T  +  1  £l 
aP”  r  3  * 


1  ill;  a<r<b,tv0 


(1) 


where  o  >  0  denotes  thermal  diffusivity.  In  addition,  the  following 
initial  and  boundary  conditions  are  specified  forT(r,t): 


T(r ,0) 

*  f(r)  ;  a<r<b,  t  =  0 

(2) 

Kil  = 

K  3r 

-h1[Tg(t)-T]  ;  r  =  a,  t  >  0 

(3) 

II 

M*- 

rc  |ro 

-h2[T-TA]  ;  r  -  b,  t  >  0 

(4) 

where  Tg( t)  is  the  propellant  gas  temperature  function  at  the  bore 
surface,  TA  is  the  ambient  temperature,  f(r)  is  the  initial  temperature 
distribution  at  the  beginning  of  each  small  time-increment,  or  the 
beginning  of  each  round,  hi  and  h2  are  respectively  gas  and  natural 
convection  coefficients,  and  K  is  the  thermal  conductivity  of  the  gun 
tube. 


The  convection  coefficients,  hx  and  h2,  at  the  bore  surface  and 
at  the  outer  surface  of  a  barrel  may  be  assumed  to  vary  at  each  small 
time-interval  during  each  firing  round.  This  justification  introduces 
a  way  to  consider  the  radiation  boundary  conditions  at  the  bore  and 
at  the  outer  surfaces  of  a  gun  tube.  The  heat  input  is  permitted  to 
vary  with  time;  in  general,  the  propellant  gas  temperature  function 
can  be  assumed  to  be  a  combination  of  various  exponential  functions. 

METHOD  OF  SOLUTION 


Since  the  boundary  condition  (3)  is  not  homogeneous,  the  method 
of  separation  of  variables  cannot  be  applied  directly.  The  Laplace 
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transform  method  is  considered;  however,  the  authors  find  that,  if 
the  Laplace  transform  technique  is  used  alone  for  solving  this  problem, 
the  inverse  transform  is  very  difficult  to  obtain.  The  combination  of 
Laplace  transform  and  separation  of  variables  will  be  adopted. 

Because  of  the  linearity  of  the  problem,  the  solution  to  this 
problem  is  readily  verified  and  may  be  written  in  the  form 


T(r,t)  =  Ta(r,t)  +  Ts(r)  +  Tb(r,t);  a  <•  r  <  b,  t  >  0  (5) 


where  Ta(r,t)  satisfies  the  following  equations: 

i2Ta  +  1  lli  =  1  lli  .  a  < 

ar7”  r  ar  a  at  * 


r  <  b,  t  >  0 


K  aF-  -hja  =  -hjg(t)  ; 

+  h2Ta  =  0  ; 

Ta(r,t)  =  0  ; 

Ts(r)  satisfies  the  following  equations: 


r  r  a,  t  >  0 

r  -  b,  t  >  0 

a-r-b,  t=0 


(6) 

(7) 

(8) 

(9) 


1 

r  dr 


0  ; 


a  <  r  <  b 


(10) 


-hJs  -  0 


r  =  a 


(ID 


dTs 

K  dT  +  Ms  =  h2T A  ;  r  -  b  (12) 
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and  Tb(rtt)  satisfies  the  following  equations: 


+  1  Ilk  =  1  |Ih.  ;  a  <  r  <  b,  t  s  0  (13) 

Dr2  r  3r  a  at  ’ 

K  ^  -hjTb  =  0  ;  r  =  a,  t  >  0  (14) 

K  +  hjb  =  0  ;  r  =  b,  t  >  0  (15) 

Tb(r,t)  =  f(r)  -  Ts(r)  ;  a^r^b,  t  =  0  (16) 


The  solution  for  Ta(r,t)  is  the  one  that  satisfies  the  homogeneous 
differential  equation  (6)  with  nonhomogeneous  boundary  condition  (7), 
homogeneous  boundary  condition  (8)  and  with  zero  initial  temperature 
distribution  (9).  The  solution  for  Ta(r,t)  will  be  obtained  by  use  of  the 
method  of  Laplace  transform.  The  problem  involving  Ts(r)  is  one  of  a 
steady-state  temperature  distribution  that  can  be  readily  solved  by  the 
usual  technique  of  solving  ordinary  differential  equations.  Once  the 
problem  regarding  Ts(r)  is  solved,  the  problem  concerning  Tk(r,t)  is  com¬ 
pletely  defined  and  is  the  one  with  homogeneous  differential  equation  (13) 
and  boundary  and  initial  conditions  (14),  (15)  and  (16).  Since  the  boundary 
conditions  involving  Tb(r,t)  are  homogeneous,  Tb(r,t)  can  be  obtained  by 
the  method  of  separation  of  variables.  When  Ta(r,t),  Ts(r),  and  Tb(r,t) 
are  obtained,  then  the  complete  solution  to  the  problem  formulated  ip 
the  previous  section  is  given  by  equation  (5). 

SOLUTION 

(A)  Solution  for  Ta(r,t) 

The  Laplace  transform  of  equation  (6)  is  given  by 


0*  +  f  Hr2-  "  92Ta  =  0  ;  a  <  r  <  b  (17) 
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where  q2  -  p/a  and  Ta  is  defined  by 


Ta 


e“pt  Ta(r,t)dt 

0 


The  complete  solution  of  the  Bessel  equation  (17)  is 


(13) 


Ta  =  Al0(qr)  +  BK^qr)  09) 


and  Kq  are  modified  Bessel  functions  of  order  zero.  The 
s  A  ana  B  are  determined  so  that  Ta  satisfies  the  Laplace 
transform  of  equations  (7)  and  (8),  namely 


where  Ip 
constant 


K 


fa 

dr 


-hJa  =  -hiTg(P)  > 


r  =  a 


(20) 


dT, 


^  dr  +Ma-0  ; 


r  -  d 


where  Tg(p)  = 


e"pt  Tg(t)dt 


(21) 


(22) 


Substituting  equation  (19)  into  equations  (20)  and  (21),  one  obtains 
A  [-hlI0(qa)  +  Kqljqa)]  -  B  DvK0(qa)  +  KqK.(qa)]  =  -hlTg(p)  (23) 


and 


A  [h2I0(qb)  +  Kql,(qb)]  +  B  [h^qb)  -  KqK,(qb)]  =  0  (24) 

Solving  A  and  B  from  equations  (23)  and  (24),  and  substituting  into 
equation  (19),  one  obtains 

T  _  -Tg(p)h,i[h2Ko(qb)  -  KqK. (qb)]i0(qr)  -  Lh^I0(qb)  +  KqI1(qb)]K0(qr)} 

j  ~  “  ’  .  _ 


A(p) 


(25) 


where  i(p)  =  [h2I0(qb)  +  Kql^qb)]  [h x K0(qa)  +  KqKx (qa) ] 


-  tVofqa)  -  Kqlj(qa)]  [h2K0(qb)  -  KqKj(qb)] 


(26) 


Suppose 


Ti  * 


-h1((h2K0(qb)  -  KqKj (qb) ] IQ(qr)  -  [h2I0(qb)  +  Kql j (qb)]KQ(qr) } 


A(0) 


Then  equation  (25)  becomes 


Ta  =  Vi 


(27) 


(23) 


T  is  now  determined  by  the  Inversion  Theorem,  i.e., 

T,  -y*  cr,]  • 


Y+i“ 


-1 

2*i 


e  pt  h1t[h2K0(qb)-KqK1(qb)3l0(qr)-[h2I0(qb)+KqI1(qb)]K0(qr))  dp 


A(p) 


(29) 


Y-i® 


The  integrand  is  a  single-valued  function  of  p  with  simple  poles  at 
p  =-a6^,  where  sn  are  the  roots  of 

[h2J0(8b)  -  KfJJjUb)]  Cn1Y0(ea)  +  KeY, (ea) ] 


ChjJ0(6a)  +  KgO. (6a)]  [h2Y0(gb)  -  KeY^gb)]  =  0 


(30) 
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Since  I0(iBnr)  =  0o(8nr),  I^iBnr)  =  U^r), 
Ko(iBnr)  =  -n/2  [iJQ(8nr)  +  Y0(enr)],  and 
K^i^r)  *  u/2  [-0,(6/)  +  iYJ6nr^ 

To  find  the  residue  at  pole  p  -  -as^ »  one  needs 


dAki  ! 

1 

dA(p) 

dp  ! 

2aX 

dA 

J 

p 

p  •  -«e* 

- 

-  x  =  l8n 

=  C  _i_ 
1“ 

f  a  [  h,!,(xa) 

+  KXI0 

, ( A.a )]  th2K0(>b)  -  KxK,(xb)] 

-  b  C-hjI0(xa)  +  Kxl^xa)]  [h/.Ub)  -  K*K0(Ab)] 


-  a  [h.I0( xb)  +  Kxl  ( xb) ]  [h  K  (xa)  +  Kx Kq ( >  a ) ] 

+  b  [h2Ij(xb)  +  KaIq( xb)]  [hlK0C^a>  -  KxKt(xa)D j  ^  (31) 

I)  '  -  I8n 

where  equation  (26)  and  the  recurrence  formulae7  have  been  used.  To 
simplify  equation  (31),  when  x  =  iBn, 

-hjl0{xa)  +  lUI^Aa)  -[h.K0(xa)  +  KaK  (xa)] 

~h2l0Ub)  *  Kxl  ( Xb)  =  hJR^UbT  nCTR  ( xb) 


-Ehl0o(6na)  +  KgnJ  (8na)] 
h2J0{,3nb)  - 


(32) 
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With  the  use  of  this  result  and  the  Wronskian  relation,  equation  (31) 
can  be  written 


I  dA(p) 
A  dA 


p( 


hj+K2^)  -  (h^K2f^) 


(h^K2S2)[-h1Jo(8na)-K0nJ1(6na)]2-(hJ+K2(32)[h2Jo(3nb)-KgnJ1(Bnb)]2 
Zap^C-hjjQtBnaJ-KSnd^Bna)]  [h2J0( Bnb)-KBn0 x ( Bnb) ] 

_ F  (Bn) _ 

2aB2[-h,J0(Bna)-K8nJ1(Bna)]  [h^Bnbl-KB^fBnb)] 


where 


F(sn)  =  (h*«*s*)  [h,J0<v)+lCi^I(iW)]MhJ+K*iiij)  Ch2 J0 ( enb ) -Ke„J x ( e„b) 32 


The  residue  at  the  pole  p  =  -aB^  or  A  =  iBn8  is 


Urn  (p+a82)  '  e  P  hi ^ Ch2K0(Ab)-KxK1(Ab)1I0(xr)-[h2I0(Ab)+KxI1(Ab)]K0(Ar) > 
•.  ift.n  Mp) 


rshi  FTbJ  6  Co^nSr)[h1J0(Bns)+XBnJ1(Bna)][h2-j0(Bnb)-i:3nJ1(Bnb)] 


wnere 


Co(en,r)  =  Ch2J0(Bnb)-KBnJ1(Bnb)]Y0(6nr)-[h2Y0(Bnb)-KBnY1(Bnb)]J0(Bnr) 


(36) 
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Thus » 


”  j-}2 

Ti(r.t)  =  "-"iJTfTFT  e  "°B"t  Coi6",r)  G(sn>  (37) 

n=i 

where 

G(  Bp }  =  [hiJ0(8na}+KBn«3i(6na)]Ch2Jo(^n^)"^n^’<(Bnb)]  (38) 

Applying  the  Inversion  Theorem  to  equation  (28),  one  obtains 

Ta(r,t)  =^'‘  [Ta]  ='x‘1  HgT,]  (39) 

Applying  the  Convolution  Theorem,8  one  notes  that 

Tgtt)  [Vg(p)3  and  Tj(r,t)  -1  [T^r.p)] 

And  finally 

t 

Ta(r,t)  =  Tg<t-0  T^ddc 

. 

o 

t 

X-—'™  -n«2  + 

=  itahj  Tg(t-s)  HTT  6  n  Co^n»r )  G(Sn) 

o  n=l  0 

^  t 

=  itahj  y~~  p^-y  CQ(en,r)  G(Bn)  Tg(t-0  e  "aSnc  dr,  (40) 

n=i  n 

o 


13 


(B)  Solution  for  Ts(r) 

The  general  solution  of  the  differential  equation  (10)  is 

T${r)  =  Alog  r  +  B  (41) 

where  A  and  B  are  constants  that  are  determined  so  that  Ts  satisfies 
the  boundary  conditions  (11)  and  (12), 

Substituting  equation  (41)  into  equations  (11)  and  (12)  and 
solving  for  A  and  B,  one  obtains 


hA  ta  _ =  h 

K  hr  -  -hih2  '°9  ?  K  4  +  t>  +  hA 109  a 
(K  ^  -hjh2  log  a)  TA 

K4’T)fhA  199  I 


(42) 


Substituting  A  and  B  into  equation  (41),  one  gets 


Mr)- 


h 

i  * 

VT 


+  -—)  r  hjh 


2  109  ? 


Chlh2 


1  °9  j  +  K  ~] ;  a  <  r  <  b 


(43) 


(C)  Solution  for  Tb(r,t) 

Once  T  (r)  has  been  obtained  from  equation  (43),  the  initial  con¬ 
dition  equation  (16)  involving  Tfcir.t)  is  completely  defined.  Since 
the  boundary  conditions  (14)  and  (15)  are  homogeneous,  the  homogeneous 
differential  equation  (13)  can  be  readily  solved  by  separation  of  vari¬ 
ables,  The  function  Tb(r,t)  =  R(r)  e(t)  is  a  solution  provided 


a  6 


(44) 
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Since  the  member  on  the  left  is  a  function  of  t  alone  and  that  on  the 
right  is  a  function  of  r  alone,  they  must  be  equal  to  a  constant,  say, 

~S*.  Hence,  one  has  equations 

rR"  +  R’  +  s2rR  =  0  (45) 

9*  +  a329  =  0  (46) 

The  equation  in  R  is  Bessel's  equation.  The  complete  solution  for 
R  can  be  written 

R(r)  =  A  Jc(Br)  +  B  Y0(S'r)  (47) 

The  solution  of  equation  (46)  can  be  written 

e(t)  =  C  e  "a82t  (48) 

where  A,  B,  and  C  are  constants  that  can  be  determined  from  equations 
(14),  (15),  and  (16).  Under  equations  (14)  and  (15),  one  finds  that 

Tb(r,t)  =  it  Dn  e  [h1Y0(B„a)  ♦  ^(8^)]  J0(snr)  (49) 

n=l 

where  a  =  6n  are  the  roots  of  the  following  equation 

k 

[h2J0(eb)  -  KsJ ( (sb) ]  [hjYgfsa)  +  KsY^sa)] 

-  [hjJotBa)  +  KeJ x ( ea) ]  [h2Y0(sb)  -  KsY^eb)]  =  0  (50) 

This  equation  can  be  noted  to  be  identical  to  equation  (30). 

At  t  =  0,  using  equations  (16)  and  (49),  one  has 

OP 

f(r)  -  T$(r)  *  21  Dn  <T0(Bn,r)  (51) 

n=l 
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where 


C0(Bn,r)  =  [hlY0(Bna)  +  KBnVM]  Jo(enr) 

-  [h  1  J0( 3na)  +  KgnJ^Bna)]  V0(enr)  (52) 

The  constants  Dn  can  be  determined  by  use  of  the  fact  that  the 
set  of  cylinder  functions 

Co(&n»n ;  n  -  1,  2)  3,  ... 

is  an  orthogonal  set  on  the  interval  a  <•  r  <  b.  Hence,  one  obtains 
b 

r  C0(en,r)  (f(r)  -  Ts(r))  dr 

Dn  =  -3 -  (53) 

b 

j  r  C"q  (en,r)  dr 
a 

In  compact  form,  equation  (49)  can  be  written 


Tb(r,t)  =  E  Dne  ’^n4  C0(pn,r)  (54) 

n=l 

where  Dn  is  defined  by  equation  (53),  and  C0(Bn,r)  is  defined  by  equa¬ 
tion  (52), 

SUMMARY 


The  analytical  solution  to  the  problem  formulated  in  a  previous 
section  is  given  by  the  sum  of  equations  (40),  (43),  and  (53).  The  gas 
and  the  natural  convection  coefficients  hj  and  h,»  respectively,  can  be 
assumed  to  be  variable  for  each  time  interval.  The  initial  temperature 
distribution  ac  the  beginning  of  each  time  interval  will  be  set  equal 
to  the  final  temperature  distribution  of  the  previous  time  interval. 

If  the  temperature  distribution  of  a  barrel  under  continuous  firing 
conditions  is  sought,  the  final  temperature  distribution  of  one  firing 
round  will  be  considered  as  the  initial  temperature  distribution  of  the 
successive  firing  round. 
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